Melting and transverse depinning of driven vortex lattices in the periodic pinning of 

Josephson junction arrays 
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We study the non-equilibrium dynamical regimes of a moving vortex lattice in the periodic pinning 
of a Josephson junction array ( JJA) for finite temperatures in the case of a fractional or submatching 
field. We obtain a phase diagram for the current driven JJA as a function of the driving current 
I and temperature T. We find that when the vortex lattice is driven by a current, the depinning 
transition at Tp{I) and the melting transition at Tm{I) become separated even for a field for which 
they coincide in equilibrium. We also distinguish between the depinning of the vortex lattice in 
the direction of the current drive, and the transverse depinning in the direction perpendicular to 
the drive. The transverse depinning corresponds to the onset of transverse resistance in a moving 
vortex lattice at a given temperature Ttr- For driving currents above the critical current we find that 
the moving vortex lattice has first a transverse depinning transition at low T, and later a melting 
transition at a higher temperature, Tm > Ttr- 



PACS numbers; 74.50+r, 74.60. Ge, 74.60. Ec 
I. INTRODUCTION 

The behavior of superconducting vortices in the pres- 
ence of periodic pinning shows very rich static and dy- 
namic phenomenalJcJ The competition between the re- 
pulsive vortex-vortex interaction and the attractive pe- 
riodic pinning potential resjilt a |i n ||Bpvel vortex struc- 
tures at low temperatures.Bfl llj'llj'tj The equilibrium 
phase transitions of these vortex structures and their 
various dynamical regimes when driven out of ^eqiiilib- 
rium are of greatJjjiterest both experimentallyuJ LrE3 and 
theoretically.ull^l~EJ Several techniques have been devel- 
oped to fabricate in superconducting samples an artifi- 
cial periodic pinningpStructure: thickness modulated s 
perconducting films, El supjercpnducting wire networks 
Josephson junction affjaySjETEl magnetic dot arrays,B sub- 
micron hole latticesJ3'B and pinning induced by Bitter 
decorationB The ground states of these systems can be ei- 
ther commensurate or incommensurate vortex structures 
depending on the vortex density (i.e the magnetic field). 
In the commensurate case, a "matching" field is defined 
when the number of vortices Ny is an integer multiple 
of the number of pinning sites Np-. Ny — nNp. A "sub- 
matching" or "fractional" field is defined when Ny is a 
rational multiple of Np-. Ny — fNp with / = p/q. One of 
the main properties of periodic pinning is that there are 
enhanced critical currents and resistance minima both 
for fractional and for matching magnetic fields, for which 
the vortex lattice is strongly pinned. „ 

In the case of Josephson junction arraysB (JJA) the 
discrete lattice structure of Josephson junctions induces 
an effective periodic pinning potential (the so-called 
"egg-carton" potential) which at low temperatures con- 
fines the_vortices at the centers of the unit cells of the 
network.! There are strong commensurability effects for 
submatching fields / = p/q, for which the vortices ar- 
range in an ordered q x q superlattice that is commen- 
surate with the underlying array of junctions. The tran- 
sition temperature Tc(/) and the critical current Idf) 



have maxima for ratiorLaL/ = p/q, which have been 
observed experimentally.H'El Moreover, as we will discuss 
here, the model that describes the physics of the JJA 
can be thought as a discrete lattice London model for 
thin film superconductors with periodic arrays of holes. 
However, this comparison can be valid only for low sub- 
matching fields since it can not describe the effects of 
interstitial vortices. 

The equilibrium phase transitions at finite tempera- 
tures of two dimensional systepis, with periodic pinning 
have been studied in the past.ll2nl3 It is possible to have a 
depinning phase transition of the commensurate ground 
states at a temperature Tp and a meltixiff Jjansition of 
the voctex lattice at a temperature TM-tlTcJ Franz and 
TeitelEil have studied this problem for the case of sub- 
matching fields. For < T < Tp there is a pinned phase 
in which the vortex lattice (VL) is pinned commensurably 
to the periodic potential and has long-range order. For 
Tp < T < Tm there is a floating VL which is depinned 
and has quasi-long-range order. For high submatching 
fields (/ > 1/30 ) both transitions coincide, Tp = Tm, 
while for low submatching fields 1/30) both transi- 
tions are different with Tp < Ta/.EJ'EJ 

The non-equilibrium dynamics of driven vortex lattices 
interacting either with random or-flecipdic pinning shows 
an interesting variety of behavior .EjE^I Many recent stud- 
ies have concentrated in the problem-of the driven VL in 
the presence of random pinning£j~Ed When there is a 
large driving current the effect of the pinning potential 
is reduced, and the nature of the fastly moving vortex 
structure has been under active discussion. The moving 
vortex phase has been proposed to be either a crystalline 
structure, a moving-glass, a moving smectic or a mov- 
ing transverse glass.EJcJ These moving phases ha\p,been 
studied both experimentallynS and numerically.Ea^El Mo- 
tivated by these results, the dynamical regimes of the 
moving VL in the presence nL.jMriodic pinning has also 
become a subject of interest.E3 At zero temperature, 
the dynamical phases of vortices driven by an external 
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current with a periodic array of pinning sites has ..been 
studied in very detail by F. Nori and coworkers .113113 A 
complex variety of regimes has been found, particularly 
for > Np where the motion of interstitial vortices 
leads to several interesting dynamical phases. 

In this paper we study the dynamical regimes of a mov- 
ing VL in the periodic pinning of a Josephson junction 
array (JJA) for finite temperatures in the case of a sub- 
matching field (/ = 1/25). We obtain a phase diagram 
as a function of the driving current / and temperature 
T, which is shown in Fig. 1. We find that when the VL is 
driven by a low current, the depinning and melting transi- 
tions can become separated even for a field for which they 
coincide in equilibrium. Moreover, we can distinguish be- 
tween the depinning of the VL in the direction of the cur- 
rent drive, and the transverse depinning in the direction 
perpendicular to the drive. This later case corresponds 
to the vanishing of the transverse critical current in a 
moving VL at a given temperature Tj^, or equivalently, 
to the vanishing of the transverse superconducting coher- 
ence. We obtain three distinct regimes at low tempera- 
tures: {i)Pinned vortex lattice: for < T < Tp{I) there is 
an ordered VL which has crystaline long-range order, su- 
perconducting coherence (i.e., a finite helicity modulus) 
and zero resistance both in the longitudinal and trans- 
verse directions, (ii) Transversely pinned vortex lattice: 
for Tp{I) < T < Ttr{I) there is a moving VL which has 
anisotropic Bragg peaks, quasi-long range order, trans- 
verse superconducting coherence and zero transverse re- 
sistivity. There is a finite transverse critical current, jlhis 
regime also has strong orientational pinning effectscil in 
the [1,0] and [0,1] lattice directions, {m) Floating vor- 
tex lattice: for Ttr{I) < T < Tm{I) there is a moving VL 
which is unpinned in both directions and it has quasi-long 
range crystalline order with a strong anisotropy. Some of 
the results we discuss here were reported,by us previously 
in a shorter and less detailed version.Ey After our work 
in Ref. |2^, further studies of thermal effects in a moving 
vortex lattice with a periodic pinning arrajif-tiave been 
reported for iV„ = iVpB and for iV„ > A^pEilS Some of 
these results are similar to ours. 

The remainder of our paper is organized as follows. In 
Sec. II we introduce the theoretical model used for the 
dynamics of the JJA. We also discuss how this model 
can be mapped to a superconducting film with a peri- 
odic array of holes. In Sec. Ill we discuss the melting 
transition near equilibrium when the JJA is driven by 
a very small current. In Sec. IV we present our results 
for the transport properties, analyzing the temperature 
dependence of the current-voltage curves and the onset 
of resistivity at Tp for different currents. In Sec. V we 
present our results for the transverse depinning transi- 
tions at different temperatures and driving currents. We 
study the transverse current- voltage characteristics both 
at T = and at finite T. We define the transverse depin- 
ning temperature Ttr from the behavior of the transverse 
resistivity. In Sec. VI we present in detail our results 
for the different dynamical regimes. For various fixed 



driving currents, we analyze the behavior of structure 
factor, longitudinal and transverse resistance and helic- 
ity modulus through the different regimes as a function 
of temperature. These results constitute the central core 
of this paper, from which the phase diagram of Fig. 1 was 
obtained. In Sec. VII we discuss the orientational pin- 
ning effects, which are observed when the direction of the 
driving current is rotated with respect to the JJA square 
lattice. Finally in Sec. VIII we present our discussion 
and conclusions. We also add an Appendix at the end, 
where we provide a detailed definition of the adequate 
periodic boundary conditions for a JJA with an external 
magnetic field and an external driving current, as well as 
the algorithm used for the numerical simulation. 

II. MODEL AND DYNAMICS 

A. RSJ dynamics 

We consider a square Josephson junction array with 
i X L superconducting nodes. The nodes are in the lattice 
sites n = {ux^ny) and their superconducting phases are 
6'(n). We study the dynamics of JJA using the resistively 
shunted junction- CESJX model for the junctions of the 
square networkE3~E3'''L3 In this case, the current flowing 
in the junction between two superconducting nodes in the 
JJA is modeled as the sum of the Josephson supercurrent 
and the normal current: 

I,in) = /osin^^n) + +^^(n,t) (1) 

ZeriN at 

where Iq is the critical current of the junction between 
the sites n and n-t-/i, (/j, = x, y), Rn is the normal state 
resistance and 

e^{n) = 0{n + li) - 0{n) - A^(n) = A^^^(n) - ^^(n) 

(2) 

is the gauge invariant phase difference with 

o p{n+fj.)a 

A,{n)^— A-dl. (3) 

^0 Jna 

The thermal noise fluctuations rj^ have correlations 

IkT 

(7y^(n,i)v(n',t')) = -^'5^,M"^n,„'^(i - t') (4) 
In the presence of an external magnetic field H we have 

X A^(n) = A.4n) - A^n + y) + Ay{n + x) - Ay{n) 
- 2^/, (5) 

/ — Ha^/^o and a is the array lattice spacing. We take 
periodic boundary conditions (p.b.c) in both directions 
in the presence of an external current I = /y in arrays 



2 



with L X L junctions.ill (See the Appendix). The vector 
potential is taken as 



A^{n,t)^Al{n)-a^{t) 



(6) 



where in the Landau gauge A°(n) = ~2tt fuy, A^[n) — 
and Q;^(t) allows for total voltage fluctuations under 
periodic boundaryi-cpaditions. In this gauge the p.b.c. 
for the phases aren3£d 

d{nx + L, Uy) = 9{nx, Uy) 

9{n,j:,ny + L) = 0{nx, Uy) - 2t: fLn^. (7) 
We also consider local conservation of current, 

• /^(n) = J2 ^p(n) - - A*) = 0. (8) 

After we obtain the following equations for the 

phasesjnaki 



where 

S^{n) = lo sin[A^0(n) - (n) - , 
and the discrete Laplacian is 



(9) 



(10) 



Alfin) = /(n + x) + /(n - x) + /(n + y) + /(n - y) 



-4/(n) 



(11) 



The Laplacian can be inverted with the square lattice 
Green's function Gn n'- 



(12) 



Since we take periodic boundary conditions (see Ap- 
pendix), the total current has to be fixed by: 



1 

l2 



^ Jo sin 6*^(11) +f]^{n,t) 



^ lo sin Oy{n)+ ijy (n, t) 



2eRN~dt 



(13) 



h day 
2eRN~dr ' 



These equations determine the dynamics of (i) -0 For 
this case we take Ix — and ly ~ I. After Eqs. (t 
we obtain the following set of dynamical equations,! 

d0{n 




dt =-Y.^r.,r.-Af.-[S^{n')+^^{n\t)\, (14) 

n' 



da 
~dt 



B. Comparison with thin film with a periodic array 
of holes 



Let us consider a superconducting thin film with a 
square array of holes, which act as pinning sites for vor- 
tices. There are Np = pinning sites separated by a 
distance a. The current density in the superconducting 
film is given by the sum of the supercurrent and the nor- 
mal current: 



J Js + Jat 

J 



— [**D*-(D*)**] + ^|- fv0-— A) (16) 



with D = V + i^A, ^'(r) = |^'(r)| exp[i6'(r)] the su- 
perconducting order parameter and a the normal state 
conductivity. These equations are valid everywhere in 
the film except in the hole regions. If the number of vor- 
tices Ny = BL^a? is much smaller than the number 
of pinning sites Np, all vortices will be centered in the 
holes in equilibrium. In this case we can assume that 
5'(n)| « Ivfol is homogeneous in the superconducting 
film. Therefore the dynamics is given by the supercon- 
ducting phase 9{v), corresponding to a London model 
in a sample with holes. After considering current con- 
servation V • J = 0, we obtain the London dynamical 
equations for the phases in this multiply connected ge- 
ometry. Since Ny <C iVp, we make the approximation of 
solving the equations in a discrete grid of spacing a. This 
means that we take as the relevant dynamical variables 
the phases 0(rn) defined in the sites which are dual to 
the pinning sites. They represent the average supercon- 
ducting phase in each superconducting square defined by 
four pinning sites. Therefore, we take the discretization 
Tn = {nxCL, Uyo) = an. [Pinning sites at centered at posi- 
tions Tp — {fix + 1/2, riy + l/2)a]. The derivatives in the 
supercurrent are discretized in a gauge-invariant way as 



i^M*(r) 



1 r 



^-(11 + ^^)- e-*2'^'4^'(")/*°*(n)l . (17) 



After doing this, we obtain an equation analogous to (|^). 
Now I^in) has to be interpreted as current density nor- 
malized by Jo = 2e?i|^'oP/TOa = $o/(87r^A^a), time nor- 
malized hy T ~ c/ (47rCTA^), and the fraction of vortices is 
/ = Ny/Np = Ba? /^Q. This leads to a set of dynamical 
equations of the same form as Eqs.( p^ , p^ . Therefore, we 
expect that for / ^ 1 the model for a JJA also gives a 
good representation of the physics of a superconducting 
film with a square array of holes (meaning that effects 
of interstitial vortices are neglected for iV„ <C Np). In 
other words, we expect that for a low density of vortices 
the specific shape of the periodic pinning potential (be- 
ing either an egg-carton or an array of holes) will not be 
physically relevant. 



where we have normalized currents by Jo, time by tj 
2eRNlo/h, and temperature by Io^o/2'KkB- 
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C. Quantities calculated and simulation parameters 

The Langevin dynamical equations (|l^,|l^) are solved 
with a second order Runge-Kutta-Helfand-Greenside al- 
gorithm with time step Ai = O.Itj. The discrete periodic 
Laplacian is inverted with a fast Fourier + tridiagonal- 
ization algorithm as in Ref. |3^. (See also the Appendix). 
We calculate the following physical quantities: 
(i) Transverse superconducting coherence: We obtain 
the helicity modulus Ta; in the direction transverse to 
the current as 



H 2v 



^sin6'^(n) 



sin 6'x(n) 



Whenever we calculate the helicity modulus along x, we 
enforce strict periodicity in by fixing oix{t) — 0. (See 
Appendix) . 

(ii) Transport: We calculate the transport response of 
the JJA from the time average of the total voltage as 



K = Mt)) = {da^{t)/dt) 



{vy{t)) = {day{t)/dt) 



(18) 



with voltages normalized by RnIq- 

(iii) Vortex structure: We obtain the vorticity at the 
plaquette n — {n^ + 1/2, Uy + 1/2) (associated to the site 
n) asLJ 



b{n) = -Af, X nint[6l,,(n)/27r] 



(19) 



with nint[2;] the nearest integer of x. We calculate the 
average vortex structure factor as 



structure factor S'(k) has the corresponding Bragg peaks 
at wavevectors G in the reciprocal space, as can be seen 
in Fig. 1(b). When the temperature is increased, the VL 
tends to disorder and above the melting temperature Tm 
a random vortex array with a liquid-like structure factor 
is obtained. Figs. 2(c) and 2(d). 

We find a single equilibrium phase transition (/ = 0) 
at Tm « 0.050 ± 0.003, which is in agreement with the 
melting temperature pbtained by Franz and Teitelcil and 
Hattel and Wheatleyy for / > 1/30. 

We now apply a very low current, / — 0.01, in or- 
der to study the near-equilibrium transport response si- 
multaneously with other quantities like structure factor 
and helicity modulus. We find a phase transition at a 
temperature Tm{I) ~ 0.046 ± 0.001, which is slightly 
lower than the equilibrium transition. In Fig. 3 (a) we see 
that there is a large jump in the resistance R = V/I 
at Tc, in good agreement adth the first-order nature of 
the equilibrium transitiono The onset of resistivity is 
a signature of a depinning transition in the direction 
of the drive. This occurs simultaneously with a melt- 
ing of the vortex lattice, corresponding to the vanishing 
of Bragg peaks, as shown in Fig. 3(b) for the two first 
reciprocal lattice vectors Gi = —(—4/25,-3/25) and 
G2 ^(-3/25,4/25). The size dependence of S(G,T) 
is shown in the inset of Fig. 3(b). In the direction of the 
current drive the helicity modulus is ill-defined since to- 
tal phase fluctuations are allowed. (See the Appendix). 
However, in the perpendicular direction to the drive the 
helicity modulus T^, can be calculated, and it is a mea- 
sure of the transverse superconducting coherence. As 
we can see in Fig. 3(c), transverse superconductivity also 
vanishes at Tm{I)- Above Tm, we find that the Tx{T) 
has large fluctuations around zero. 



IV. TRANSPORT PROPERTIES 



j2 X! exp(ik • ii) 



(20) 



III. TRANSITION NEAR EQUILIBRIUM 

We study JJA with a low magnetic fleld corresponding 
to / = 1/25, for different system sizes oi L x L junctions, 
with L = 50, 100, 150. Most of the resuhs are for L = 
100, except when it is explicitly specified, and for Nt = 
10^ iterations after a transient of Nt/2 iterations. 

The ground state vortex configuration for / = 1 /25 is a 
tilted square-like vortex lattice (VL)J13 see Fig. 2(a). We 
find that this state is stable for low currents and low tem- 
peratures (in fact, the structure of Fig. 2(a) corresponds 
to / = 0.01 and T = 0.01). The lattice is oriented in the 
[4a, 3a] direction and commensurated with the underly- 
ing periodic pinning potential of the square JJA. The 



Let us now study the transport properties for larger 
currents. We calculate the current-voltage (IV) charac- 
teristics for different temperatures as well as the dc re- 
sistance R — V/I as a function of temperature for finite 
currents. 

The zero temperature IV curve has a critical current, 
/c(0) = 0.114 ± 0.002, see Fig. 4, which corresponds to 
the single vortex depinning current in square JJAJ3 with 
the typical square root depence at the onset. Similar be- 
havior has been reported for zero temperature IV curves 
for low values of /.EHu Above /c(0) there is an almost 
linear increase of voltage, corresponding to a "flux flow" 
regime, where there is a fastly moving VL. The structure 
factor of the T = moving VL is the same as the corre- 
sponding one of the pinned VL [Fig. 2(a)]. The presence 
of periodic boundary conditions in our case prevents the 
occurrence of random or chaotic vortex motion near the 
critical current, as reported in early simulations with free 
boundary conditionsHaS In what follows we will restrict 
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our analysis for currents / < 0.4, where the collective 
behavior of the VL is the dominant physics. (At 1^1 
there is a sharp increase of voltage when all the junctions 
become normal, and V ^ RnI for / ^ 1). 

The IV curves for finite temperatures are shown in 
Fig. 5. For temperatures below Tm there is a nonlinear 
sharp rise in voltage which defines the apparent critical 
current Ic{T). For example, we can obtain this Ic{T) 
with a voltage criterion, which we choose asV < j^TAt ~ 
10"^. In this case, we find that Ic{T) decreases with T, 
vanishing at Tm- It is interesting to point out that all the 
IV curves for different temperatures have a crossing point 
at /* = 0.165, see Fig. 5(a). A crossing in the IVs has alsa 
been reported in experiments in amorphous thin films.cll 
For temperatures T > Tm the IV curves tend to linear 
resistivity for low currents. This is shown in Fig. 5(b) in 
a log-log plot of R{I) — y// vs /, where we see that R{I) 
tends to a low current finite value for T > Tm, while it 
has a strong nonlinear decrease for T < Tm- 

Let us now study the dc resistance R = V/I as a func- 
tion of temperature for a given applied current in the 
y-direction (Fig. 6). We start with the perfectly ordered 
VL as an initial condition at T = and then we slowly 
increase the temperature, keeping / constant. For cur- 
rents below the T = critical current, / < /c(0), the dc 
resistance is negligibly small at low T, and it has a steep 
increase at a depinning temperature Tp(I), corresponding 
to the onset of vortex motion. The depinning tempera- 
ture decreases for increasing currents, and the values of 
Tp(I) are coincident with the apparent critical currents 
Ic{T) obtained from the IV curves. For currents higher 
than /c(0), there is always a large and finite voltage for 
any temperature. If /c(0) </</*, the R{I) increases 
slightly with T tending to a constant value for large T, 
while for I > I* the R{I) decreases with T. 

V. TRANSVERSE DEPINNING 

What is its response to a small current in the transverse 
direction when the driven vortex lattice is moving? Is the 
vortex lattice still pinned in the transverse direction? Is 
there a transverse critical current for a moving VL? The 
idea of a transverse depinning current was introduced by 
Giamarchi and Le Doussal in Ref. for moving vortex 
systems with random pinning at zero temperature. The 
possibility of such a critical current was later questioned 
by Balents, Marchetti and Radzihovsky^H where it was 
shown that this is not true for any finite temperature 
in random pinning; however a strong nonlinear increase 
of the transverse voltage was predicted at an "effective" 
transverse critical current. In the case oi periodic pinning 
it is more clear that a transverse critical current will exist 
at T = since it is a commensurability effect. This has 
been fmmd in the T = simulation work of Reichhardt 
et al.x3 It is also possible that this transverse critical 
current will still be non-zero at T ^ in periodic pinning. 



In fact, we have found in our previous workE3 that there 
is a thermal transverse depinning in a periodic system, 
and we will now analyze this behavior in detail. 

First, a high longitudinal current / > /c(0) is applied 
at zero temperature. Then, a current Itr is applied in the 
transverse direction. In this way, a transverse current- 
voltage characteristics can be obtained for each /. This is 
shown in Fig. 7(a) for / = 0.16. We clearly see that there 
is a finite transverse critical current Ic,tr ~ 0.12, which is 
of the order of the single vortex pinning barrier. We also 
show in Fig. 7(a) how the longitudinal voltage V changes 
when the transverse current is varied. For Itr < Ic,tr the 
longitudinal voltage V is almost constant. At Ic.tr there 
is a fast decay of V. When Itr — I we also have V = Vtr 
as expected for a drive at a degree of 7r/4. Later, for 
Itr > I the vortex lattice becomes pinned in the other 
direction, since now the directions of "longitudinal" and 
"transverse" current are interchanged. 

Another possible measurement is to study thermal 
transverse depinning. In this case, we start with a lon- 
gitudinal current /, then a small transverse current is 
applied, Itr ^ Ic.tr, and the temperature is slowly in- 
creased. In this way, we can measure a transverse re- 
sistance Rtr = Vtr /Itr- In Fig. 7(b) we plot this result 
for / = 0.16 and Itr = 0.01. We find that for finite low 
temperatures Rtr is negligibly small within our numerical 
accuracy and it has a clear onset at a transverse critical 
temperature Ttr- 

Let us now see how these results depend on the longi- 
tudinal current / and temperature T. We have calculated 
the transverse IV (Tr-IV) curves for different I and T. In 
Fig. 8(a) we show the Tr-IV curves for / = 0.06 (low cur- 
rent regime) and in Fig. 8(b) for / = 0.16 (high current 
regime). In both cases, there is a clear change of behavior 
in the Tr-IV curves when going through a characteristic 
transverse critical temperature Ttr{I). For low tempera- 
tures there is a transverse critical current which tends to 
vanish when T aproachs Ttr{I) from below. In contrast, 
for T > Ttr there is a linear resistivity behaviour. 

The transverse resistivity Rtr = Vtr I Itr as a function 
of temperature was calculated for different longitudinal 
currents (Fig. 9). In all cases there is an onset of trans- 
verse response at a given temperature Ttr{I). At low cur- 
rents, I < /c(0), the transverse depinning temperatures 
are almost constant, Ttr ^ 0.02 tending to increase slowly 
with /, see Fig. 9(a). On the other hand, for / > /c(0) the 
transverse depinning temperatures increase clearly with 
/, see Fig.9(b). 

VI. NON-EQUILIBRIUM REGIMES 

We will now study the different non-equilibrium 
regimes of vortex driven lattices and characterize their 
possible dynamical transitions. The approach we will 
follow in this Section is to have a fixed current applied 
in the system and vary the temperature. In this way, we 
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look for the possible transitions as a function of T in a 
similar way as was done near equilibrium in Sec. III. 

A few similar studies were done previously in related 
systems. In Rcf. ^ the melting transition of a moving 
vortex lattice in the three dimensional XY model was 
studied in this way. In this work a first order transi- 
tion was found as a function of temperature in a strongly 
driven vortex lattice. In Ref. ^ a current driven two 
dimensional JJA at zero field was studied. The possibil- 
ity of a transition as a function of temperature for finite 
currents below the critical current was analyzed in this 
case. 

In the following, we will separate our study in three 
ranges of current: (i) low currents, / < 0.04, (ii) inter- 
mediate currents 0.04 < I < /c(0), and high currents 
/ > /c(0). 

A. Low currents 

We show the results for a low current in Fig. 10 for 
/ — 0.03. The longitudinal dc response, V/I, is negligi- 
bly small at low T and later it has a sharp increase of two 
orders of magnitude: this defines the depinning temper- 
ature, Tp [Fig. 10(a)]. At higher temperatures, T > Tp, 
the resistance is weakly T-dependent. 

Below Tp the vortex lattice is pinned and its structure 
is similar to the ground state: a vortex lattice commen- 
surate with the underlying square array and tilted in the 
[4a, 3a] direction. Above Tp the vortex lattice is moving 
and it has an anisotropic structural order. If we ana- 
lyze the structure factor in two different reciprocal lat- 
tice directions S{Gi) and S{G2), we can see this clearly 
(Fig. 10(b)). Below Tp the VL structure is isotropic and 
S{Gi) = S{G2)- Right above Tp the height of the peaks 
decreases with temperature, and the structure of the de- 
pinned VL is clearly anisotropic, S{Gi) ^ S{G2)- Fi- 
nally, at a melting temperature Tm the peaks vanish, 
and the vortex lattice melts into a vortex liquid. In 
Fig. 11 we compare the behavior of the Bragg peaks for 
two system sizes L = 50,100. We see that the S{Gi,2) 
are size independen±_.below Tp as it should be expected 
for a pinned phase.Ej For large temperatures T > Tm, 
in the liquid phase, the value of S'(G'i_a)-|is strongly size 
dependent, since it should go as On the other 

hand, for Tp < T < Tm we find that the intensity of the 
Bragg peaks is weakly dependent on system size. The 
clear change of behavior of the size dependence gives a 
good criterion to determine Tm (Fig. 11). 

The helicity modulus in the direction perpendicular to 
the current, T^;, decreases very slowly for T < Tp. Above 
Tp, Tj; has a faster decay with important fluctuations and 
tends to vanish at Tm- For T > Tm, '^x oscillates around 
zero. 

Therefore, for a small finite current the depinning and 
melting transitions become separated with Tp < Tm- 



B. Intermediate currents 

At intermediate currents, 0.04 < I < /c(0), a new tran- 
sition appears: the transverse depinning of the moving 
vortex lattice. As discussed in Sec.V, one can measure 
transverse depinning by applying a small transverse cur- 
rent while the VL is driven with a fixed longitudinal cur- 
rent. This is shown in Fig. 12(a) for the case of / = 0.06 
and a small transverse current, Itr — 0.01. We see that 
there is an onset of transverse voltage at Ttr — 0.02. We 
can also see that this transverse depinning temperature 
is above the depinning temperature Tp for longitudinal 
resistance [Fig. 12(b), Tp = 0.013], and below the melt- 
ing temperature Tm for the vanishing of the Bragg peaks 
[Fig.l2(c), Tm = 0.037]. Therefore, this transition oc- 
curs at an intermediate temperature between the depin- 
ning and the melting transitions, Tp < Ttr < Tm- We 
also show in Fig. 12(d) that the helicity modulus begins 
to fall down slowly at Ttr, while for Ttr < T < Tm it has 
strong fluctuations, being difficult to interpret its behav- 
ior in this case. 

It is interesting to study in more detail the behavior 
of the structure factor through all these transitions. The 
temperature dependence of the intensity of the Bragg 
peaks is shown in Fig. 12(c). In Fig. 13 we show examples 
of the structure factor S'(k) at temperatures in the dif- 
ferent regimes. In the pinned phase, the <S'(k) is nearly 
the same as in the ground state with delta-like Bragg 
peaks, see Fig. 13(a). For the moving VL, we can see that 
there is less anisotropy in the transversely pinned regime 
Tp < T < Ttr [Fig. 13(b)] than in the floating regime 
Ttr <T < Tm [Fig.l3(c-d)]. Also in the inset of Fig. 12(c) 
one can see that the intensity of the Bragg peaks has a 
greater dependence with size for Ttr < T < Tm when 
compared with the Tp < T < Ttr regime. Moreover, near 
Tm the VL structure becomes strongly anisotropic, with 
the peak at Gi much larger than the peak at G2, see 
Fig. 13(d). 

The anisotropy of the Bragg peaks of the moving VL 
(in the regimes at Tp < T < Tm) has two characteristics: 
(i) the width of the peaks increases with T in the direction 
of the applied current (the direction perpendicular to the 
vortex motion) , and (ii) the height of the peaks decreases 
in the direction of vortex motion. This can be observed 
in the sequence of structure factors shown in Fig.l3(b- 
d) . The details of the broadening of the Bragg peaks are 
shown in Fig. 14(a) for the transversely pinned regime 
(Tp < T = 0.018 < Ttr) and in Fig.l4(b) for the floating 
solid regime {Ttr <T = 0.023 < Tm)- 

In Fig. 15 we show the temperature dependence of 
S'(Gi) and S{G2) for different currents. In all the cases 
the three different regimes in temperatures are observed: 
(i) pinned regime for < T < Tp with isotropic and size 
independent Bragg peaks, (ii) transversely pinned regime 
for Tp < T < Ttr with weakly anisotropic Bragg peaks, 
and (iii) the floating solid regime for Ttr < T < Tm 
with strong anisotropy and size dependence in the Bragg 
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peaks. 



VII. ORIENTATIONAL PINNING EFFECTS 



C. High currents 

In the case of currents larger than /c(0), the vortex 
lattice is already depinncd at T = 0. As we have seen 
in Sec. V, this zero-temperature moving VL has a finite 
transverse critical current, and therefore it is pinned in 
the transverse direction. When we slowly increase tem- 
perature from this state, we find that the transverse 
resistive response is negligible for finite low tempera- 
tures. At a temperature Ttr there is a jump to a fi- 
nite transverse resistance Rtr = Vtr/Itr- For example, 
this is shown for / = 0.16 with a small transverse cur- 
rent, Ix = 0.01 in Fig. 16(a). The vortex lattice has 
an anisotropic structural order for all temperatures, i.e, 
S'(Gi) 7^ S'(G2), and the height of the Bragg peaks van- 
ishes at Tif, Fig. 16(b). In the inset we show S'(G2) for 
different sizes, L = 50, 100, and we see that Tm is size 
independent. Similar behaviour is found for S'(Gi). The 
transverse helicity modulus (T^) is almost constant for 
low temperatures and starts to decrease at Ttr presenting 
strong fluctuations for T > Ttr , Figl6(c). 

The analysis of the heights of the peaks in S'(k) as a 
function of system size L is a good indicator of the trans- 
lational correlations in the system. This_dependence is 
well known for two dimensional lattices. til For a pinned 
solid S{G) ~ 1, for a floating solid, S'(G) ~ L-'^oiLT) 
with < riQ{I,T) < 2, being this dependence a signa- 
ture of quasi- long range order, and for a normal liquid 
S{G) ~ L~^. To assure the existence of algebraic trans- 
lational correlations, we have done this scaling study for 
currents 0.02 < I < 0.2 and different temperatures for 
system sizes of L = 50, 100, 150. For the cases corre- 
sponding to the pinned regime (T <Tp) we found Tyc ~ 0, 
as expected. In the transversely pinned regime, we show 
a case in Fig. 17 (a), we find a power law fitting with very 
small values of 770(1, T). In the floating regime, we flnd 
larger values of j]g, we show a case in Fig. 17(b). In all 
the cases we have obtained that < rja {I, T) < 2 for 
Tp < T < Tm- Therefore, this flnite size analysis shows 
the existence of quasi-long range order in the moving VL. 
Also, we find that ryc^ > r]G^ for all currents and tem- 
peratures. The power-law exponent rjc can be studied at 
a constant current, as a function of temperature. This 
is shown in Fig. 18, for / = 0.16 and different reciprocal 
lattice vectors. The exponent ry is finite for the complete 
temperature range. For T < Ttr it has a small value 
rjG ~ 0.01. It has a fast increase near Ttr, where there is 
also a clear difference between r]Gi and 77^2 . Finally, it 
reaches a value of 7; = 2 near T — Tm ■ 



A very interesting characterization of the different 
regimes can be obtained by, studying the effects of vary- 
ing current direction.OE3c3 We apply a current / at an 
angle (j) with respect to the [10] lattice direction, 

Ix = I cos (j) 

Iy=Isin(j). (21) 

We study the voltage response when varying the orienta- 
tion (p of the drive while keeping fixed the amplitude / of 
the current. In the parametric curves of Vy{4>) vs. V^((/>) 
can we analyze the breaking of rotational symmetry in 
the different regimes of / and T. In the case of rota- 
tional symmetry this kind of plot should give a perfect 
circle. However, the square symmetry of the Josephson 
lattice will show up in the shape of the curves. In Ref. |2^ 
the motion of single vortices was studied. In this case, 
the "diagonal" [11] direction is unstable against small 
changes in the angle (j), while the [10] and [01] direc^. 
tions are the preferred directions for vortex motion.Efl 
This shows as "horns" in parametric vs. T4 plots, 
which are finite segments of points lying in the x or the 
y axis. This implies the existance of transverse pinning 
in these directions (thus, it corresponds to orientational 
pinning). Here, we perform the same analysis for the 
different regimes of the moving VL. 

In Fig. 19 we plot the voltages Vy and Vx when varying 
the orientational angle (f) for different current amplitudes 
/ and temperatures T. In Fig. 19 (a) we have / — 0.06 and 
T = 0.015, corresponding to the regime of a transversely 
pinned lattice. In this case most of the points are ly- 
ing either on the axis 14 = or on the axis Vy — 0, 
indicating strong orientational pinning in the symme- 
try lattice directions [10] and [01]. When increasing T 
the orientational pinning decreases and the length of the 
"horns" in the x and y axis decreases. Fig. 19(b) shows 
results for / = 0.06 and T = 0.03, which correspond to 
T > Ttr when there is a finite transverse resistance. In 
this case the horns have disappeared and orientational 
pinning is lost. However, the breaking of rotational sym- 
metry is still present in the star-shaped curves. Also in 
the high current regime. Fig. 19(c), for a low temperature 
T — 0.015 (T < Ttr) we see that there is orientational 
pinning with the presence of horns, which again dissap- 
pear for T > Ttr as it is shown in Fig. 19(d) at T = 0.035. 
Finally for T ^ Tm, deep inside in the liquid phase, the 
stars tend to the circular shape of rotational invariance. 

VIII. DISCUSSION 



With all the information of the previous sections we 
obtain the current-temperature phase diagram which was 
advanced in Fig.l. For finite currents we have been able 
to identify three different regimes, a pinned VL for T < 
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Tp{I), a transversely pinned VL for Tp{I) < T < Ttr{I), 
and a floating VL for Ttr{T) < T < Tm{I)- It is, how- 
ever, difficult to define if the temperatures Tp,Ttr,TM 
correspond to either phase transitions or to dynamical 
crossovers. The comparisons we have made of the be- 
havior of voltages (longitudinal and transversal), struc- 
ture factor and helicity modulus show that something is 
happening at these temperatures. Also, the comparison 
for different system sizes of the behavior of S{G) suggest 
transitions for Tp, Ttr, Tm- The transverse helicity modu- 
lus has strong fluctuations for Ttr < T < Tm ■ These fluc- 
tuations are not reduced when increasing the simulation 
time in a factor of 10. This could mean that actually the 
range of temperatures of Ttr < T < Tm corresponds to 
a long crossover region towards a liquid state. Also, one 
may question if the use of the helicity modulus for these 
far from equilibrium states is correct, since T^; has been 
defined from the response of the equilibrium free energy 
to a twist in the boundary condition. Most likely, the 
large fluctuations in Ttr < T < Tm are due to the very 
unstable and history-dependent steady states we flnd in 
this regime. For example the steady states obtained when 
increasing temperature from T = at flxed current differ 
from the steady states obtained when increasing current 
at fixed temperature in this regime. They differ in the 
degree of anisotropy and intensity of the Bragg peaks in 
the structure factor. We think that this reflects the fact 
that the VL is unpinned in all directions and the orien- 
tation of the moving state will depend on the history for 
this regime. 

In any case, we can clearly distinguish different dynam- 
ical states which deflne the regimes shown in the phase 
diagram of Fig.l. From the experimental point of view 
these transitions (or crossovers) are possible to measure. 
The depinning temperature Tp can be obtained from re- 
sistance measurements at finite currents. The transverse 
current-voltage characteristics and the transverse resis- 
tivity can be measured for different longitudinal currents 
and temperatures, and therefore Ttr could be obtained. 
The helicijtv modulus can be measured with the two-coil 
techniqucE It could be very interesting to see the re- 
sults of such measurements at finite currents. The struc- 
ture factor and melting transition can not be measured 
directly. However, in the presence of an external rf cur- 
rent, the dissappeatance of Shapiro steps can indicate the 
melting transition,EZI since they as^ensitive to the trasla- 
tional order of the vortex lattice.cS Therefore, we expect 
that the results obtained here could motivate new exper- 
iments for fractional or submatching fields in Joscphson 
junction arrays as well as in other superconductors with 
periodic pinning. 
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APPENDIX A: PERIODIC BOUNDARY 
CONDITIONS 

1. Phases 

We want to obtain the periodic boundary condition 
(PBC) for superconducting phases 9{nx,ny). In general 
we can write the PBC as: 

6'(n + Lx) = 9{n)+u^{n) 

e{n + Ly) = 9{n)+Uy{n) , (Al) 

where Lx = {Lx,0) and Ly = (0, L^). Taking into ac- 
count that all variables should be independent of the or- 
der of subsequent global translations and that the phases 
are defined except for an addition of 27rZ, we are lead to 
the consistency condition 

Ux{n) + Uy{n + Lx) = Ua;(n + Ly) + Uy{n) + 2iil , (A2) 

Therefore, in order to specify the periodic boundary con- 
ditions we have to give the functions Ux,Uy, which will 
depend on the gauge for the vector potential A. 

The periodic boundary condition for the phases can 
be deduced by requiring that all physical quantities are 
invariant after a translation in the lattice size. The su- 
percurrents Sp{n) = sin[Ap6'(n) — ^^(n)] should satisfy: 

^^(n + Lx) = ^^(n) 

5^(n + Ly) = ^^(n) . (A3) 

This implies that the gauge invariant phase difference 
6'p(n) = Ap6'(n) — Ap{n) should satisfy 

0^(n + Lx) = 0^(n) + 27r; 

6p{n + l.y) = 6p{n) + 2^l' , (A4) 

with I, V any integer. This condition leads to 

A^'U:^(n) = 2tiI + Ap{n + Lx) - Ap{n) 

A^uy(n) = 2ttI' + A^(n + Ly) - ^^(n) . (A5) 

We can choose the solution with I = V = {]. In the 
Landau gauge. Ax is a linear function of Uy and Ay is a 
linear function of rix ■ Taking the origin such that A^ (n = 
0) = 0, we obtain for the Landau gauge: 

Ua,(n) = Ux{Q) + Ay (Lx)7ij, 

Uy{n) = Uy{Q) + Ax{l^y)nx ■ (A6) 

The consistency condition (A2) requires 

Ay{L^)Ly- Ax{Ly)Lx=2TTl . (A7) 
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The term in the left side is equal to the total flux 
2Tr fL^Ly, therefore (A7) is equivalent to flux quantiza- 
tion, giving fLxLy = Ny, with Ny = I the number of 
vortices. 

If we take Ux{0) = Uy{0) = 0, we obtain for the PBC 

9{nx + Lx, Uy) = 9{nx,ny) + Ay (Lx)% 

0{nx, Uy + Ly) = e{nx, ny) + Ax{l,y)nx. (A8) 

A particular choice can be the gauge with Ax(n) = 
—2'jTfny, Ay{n) = 0, which leads to Eq. (7). 



2. External currents and electric fields 

In the presence of external currents or voltages the pe- 
riodic boundary conditions have to be reconsidered. In 
this case it is possible to have E • dl 7^ in a path 
that encloses all the sample either in the x or the y di- 
rection. Therefore, in a closed path we have to consider 
the Faraday's law 



The average electric field in the fj, direction is: 



dt 



E • dl = - 



c dt 



The two-dimensional sample with PBC can be thought 
as the surface of a torus in three-dimensions. The closed 
paths we are considering are the two pathes that encircle 
the torus. The electric field is not a gradient of a poten- 
tial, it is now given by E = —W — One possible 
solution is to consider a vector potential 

A(r,t) = Ao(r)-a(t), 

for which 

Hext = V X Ao 

Ida 

Etotal = ~~7r 
c dt 



where we have used the fact that A^^^^ = (which 
is the discrete equivalent of f W • dl = 0). The current 
in the link n, is, in normalized units,: 

with S'^(n) = <S'/n(n) +ry^(n, f). Therefore the average 
current in the fi direction is 

There are two cases to consider: (i) external current 
source: the external current is given and the total voltage 
fluctuates and (ii) external voltage source: the external 
voltage is given and the total current fluctuates. 

(i) External current source: the average current in /x 
direction is flxed by the external current: If?' = If^^. 
The average electric field is a fiuctuating quantity given 
by: 



Ef 



rW = ^=r-I^EW (AlO) 



In this case, the Q;^(i) is a dynamical variable, and its 
time evolution is given by Eq.(AlO). 

(ii) External voltage source: the average electric field 
in the /i direction is fixed by the external electric field: 
= E'^*'. Therefore, now the is given by the ex- 
ternal source: 



a^{t) = I E;-'dt, 



as: 



In our case, we take the adimensional vector potential ^nd a^(t) = E^^H, if E"^^ is time-independent. The 

average current is now a fluctuating quantity given by: 



A^in,t)=Al{n)-a^{t), (A9) 

with A^{n) in the Landau gauge {A^{n) = — 27r/%, 
Ay{n) = 0). Therefore the gauge invariant phase is: 

e^{n, t) = A^e{n, t) - A° (n) + a^{t). 

Then acts as a global time-dependent phase in the n 
direction. 

In the normalized units used in this paper, the electric 
field in the link defined by the junction n, n is 



Let us see how the PBC are affected by a change of 
gauge. The gauge transformations are the following: 

6»(n) ^ 61' (n) = e{n) + (l){n) 
A^{n) ^ A'^{n) = A^{n) + A^0(n) 

#(n) 



V{n) V'{n) = V{n) - 



dt 



dA^jn) 
dt ' 



where the electrostatic potential is V{n) 
Therefore we have 

' dt dt 



rf6>(n) 
dt ■ 



An interesting choice is: 

^(n) = a • n = ^ a^n^ 

In this gauge we have 6*^(11) = An9{n) — A°(n) and the 
PBC for the phases is: 
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6'(n + Lx) = 6'(n) + Ay(L^)ny + a^L^ 
0{n + Ly) = e{n) + Ax(Ly)nx + ayLy, 



and for the voltages: 
Vin + -L^)^V{n) 



da 



dt 



V{n) 



(All) 



(A12) 



The equations of motion in this gauge are: 



dejn) 
dt 

da^ 
~dF 



da 



jCKt _ 



LxLy 



(A13) 



The periodic boundary conditions with a fixed external 
current using Eq. (AlO) was used previously in Ref. ^ 
for / = (it was called a "fluctuating twist boundary 
condition") and in Ref. |3| for / 7^ 0. Also, the periodic 
boundary conditions in the gauge of Eqs. (A11-A12) were 
used in Ref. ^ for a time dependent d-wave Ginzburg- 
Landau model. 



3. Helicity modulus 

The helicitity modulus expresses the "rigidity" of 
the system with respect to an applied "twist" in the pe- 
riodic boundary conditions. The twist A:^ is defined as a 
phase change of L^k^^ between the two opposite bound- 
aries which are connected through the PBC in the n di- 
rection 



9{n + -L^) = din) 



The helicity modulus is obtained from the free energy 
F{T, k^) as: 



T — 



1 d^FiT,k^) 



dki 



(A14) 



It is clear from Eq. (All) that k^ — a^. Then, in or- 
der to evaluate the helicity modulus, a^ must be set to 
zero. This means that the helicity modulus can not be 
calculated in the direction in which there is an applied 
current, since it gives a fluctuating twist fc^(t). 



/(n-|-Lx) = /(n) 
/(n + Ly) = /(n). 



(A16) 



The linear system of L^Ly equations of ( A_15) is singular. 
Physically, the reason is that in Eq.(^), /(n) = dd{n)/dt 
corresponds to a voltage, which is defined except for a 
constant. We choose the voltage reference such that it 
has zero mean: 







(A17) 



other choices are also possible (like, for example, fixing 
/(no) = at a given site no). 

The method we use to invert the Eq. ( A15| ) is based 
on the Fourier Accelerated and Cyclic Reduction (FACR) 
algorithm.c3 In this case, we take first the discrete Fourier 
transform in the x direction: 



d{kx 1 '^y 



d{nx,ny)e 



(A18) 



which with a fast Fourier algorithm takes a computation 
time of order log Lx ■ This leads to the following equa- 
tion: 

^Kf{kx,ny) - f{kx,ny - 1) - f{kx,ny -|- 1) = d{kx,ny) , 

(A19) 

with =4 — 2 cos and with boundary condition 

f{kx,Ly + l) = /(fcj:, 1). This is a cyclic tridiagonal equa- 
tion which can be solved with a simple LU deconaposition 
algorithm in a computation time of order Lyr3 In this 
way the f{kx,ny) is obtained from (A19). Finally we 
take the inverse Fourier transform to obtain 



f{nx,ny) = -^^Rkx, 



ny)e~ 



(A20) 



This algorithm takes a computation time which is of or- 
der LxLy{logLx + A) with the constant A ^ 1. This 
is faster than the two dimensional Fourier tcansform 
method used by Eikmans and van HinibergenJ^a which 
takes a computation time of order L^; Ly (log La; + log Ly). 



4. Algorithm 

The core of the numerical calculation is to invert the 
Eq. (9). This means to solve a discrete Poisson equation 
of the form, 

A2/(n) = d{n) (A15) 
with the periodic boundary conditions 
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FIG. 1. I -T Phase diagram for / = 1/25. Tm(I) line is 
obtained from S{G) vs. T curves (•). Tp{I) line is obtained 
from IV curves (★), from S{G) vs. T curves (•) and from (Vy) 
vs. T curves (A). Ttr{I) curve is obtained from (Vtr) vs. T 
curves (♦). 
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FIG. 2. Vortex configuration 

for a low current I = 0.01 <C lc{0) and a low temperature, 
T — 0.01 < Tm: (a) tilted square vortex lattice, oriented in 
the [4a, 3a] direction and (b) corresponding structure factor 
S'(k). At high temperature, T = 0.05 > Tm- (c) disordered 
vortex array, (d) vortex liquid-like structure factor. 
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FIG. 3. For low current / < /^(O), / = 0.01: (a) 
< Vy > /ly vs. T, (b) 5(Gi) {o) and ^(Ga) (★) vs. T, 
inset: size effect in 5(Gi), (c) Tj, vs. T. 
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FIG. 4. Zero temperature IV curve. There is a critical 
current 4(0) = 0.114 ± 0.002. 
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FIG. 5. IV curves for different temperatures, (a) Lin- 
car scale, erossiug poiut at I* = 0.165. (b) Log-log plot of 
R = V/I vs. / for the same temperatures. 
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FIG. 7. Transverse depinnig temperature determiriatioris: 
(a) Tr — IV curve and < Vy > vs Itr/Iy at T = and high 
longitudinal current, ly = 0.16. (b) With small transverse 
current applied, Ix = 0.01, and ly = 0.16, < Vtr/hr > vs T. 
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FIG. 6. < Vy > /ly vs temperature curves for different dc 
currents (7^). 
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FIG. 8. Transverse IV curves for different temperatures: 
(a) low current ly < /c(0), ly = 0.06, (b) high current 
ly > /c(0), ly = 0.16. 
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FIG. 9. Vtr/hr vs T curves for different dc currents ap- 
plied in the y direction and small transverse current applied, 
Ix = 0.01: (a) low current ly < /c(0), (b) high current 

ly > /c(0). 
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FIG. 10. For low current I < 7c (0), I = 0.03: (a) 
< Vy > /ly vs. T, (b) 5(Gi) («) and ^(Ga) (*) vs. T, 
(c) vs. T. 
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FIG. 11. For low current I < 7c(0), I = 0.03: size effect in 
^(Gi). 



14 



0.0100 



0.0010 r 




0.01 0.02 0.03 0.04 0.05 

Temperature 



FIG. 12. For / < /^(O), / = 0.06: (a) < Vtr > /hr vs. T, 
(b) <Vy> /ly vs. T, (c) S'(Gi) (o) and ^(Ga) (★) vs. T, 
inset: size effect in S'(Gi), (d) T^, vs. T. 




FIG. 13. Intensity plot of the structure factor S{G) at 
I < 7c(0), / = 0.06 for different temperatures: (a) T < Tp, 
T = 0.0005, (b) Tp<T < Ttr, T = 0.015, (c) Ttr<T< Tm, 
T = 0.025, {d)T <Tm,T = 0.034. 
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FIG. 14. Details on the structure factor 
peak, Gi = ^(-4/25,-3/25), at / < 7^(0), I = 0.06, for 
different temperatures: (a) Tp < T < Ttr, T = 0.018, (b) 
Ttr <T <Tm,T = 0.023. 
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FIG. 16. For high current / > 7c(0), / = 0.16: (a) 
< Vtr > /Itr vs. T, (b) ^(Gi) (o) and S{G2) (*) vs. T, 
inset: size effect in S{G2), (c) vs. T. 
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FIG. 15. S(Gi) (o) and 5(G2) {*) vs. T for different dc 
current applied in y direction: (a) ly = 0.04, (b) ly = 0.05, 
(c) 7^ = 0.06. 
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FIG. 17. Finite size analysis and power law fit of 
S{G) ~ L-'i<^{I,T) for a moving vortex lattice at 7 > 7^(0), 
I = 0.16. We obtain: (a) for T <Ttr,T = 0.02, ?7Gi = 0.023 
(★) and J7G2 = 0.013, (A), (b) for Ttr <T <Tm,T = 0.035, 
rjGi = 0.471 (★) and tjgj = 0.34 (A). 
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FIG. 18. At 7 = 0.16, r){I, G) vs T, for Gi and G2. 
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